Realization of Geometric Landau-Zener-Stiickelberg Interferometry 
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We report the first experimental realization of the geometric Landau-Zener-Stiickelberg (LZS) 
interferometry proposed by [Phys. Rev. Lett. 107, 207002 (2011)] in a single trapped ion system. 
Different from a conventional LZS interferometer, the interference fringes of our geometric interfer- 
ometer originate solely from geometric phase. We also observe the robustness of the interference 
contrast against noise or fluctuation in the experimental parameters. Our scheme can be applied to 
other complex systems subject to relatively large errors in system control. 
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Any quantum operation, performed for quantum infor- 
mation processing or the control of quantum dynamics 
in complex systems, is subject to two primary kinds of 
errors: (i) random phase errors arising from environmen- 
tal coupling or (ii) operational errors originating from 
imperfections or fluctuations in experimental control [1|. 
Recently there have been extensive theoretical and ex- 
perimental studies in reducing random phase errors by 
the use of dynamical decoupling ]^\^ alongside quantum 
error corrections [a, Iq]. In many quantum systems, how- 
ever, operational errors may play a more significant role 
in the fidelty of the relevant quantum operations. Com- 
posite pulses or adiabatic manipulation with geometric 
phase have been extensively studied to reach error rates 
below the fault-tolerant level with reasonable limitations 
of control in feasible physical systems [7|, |8| . 

Generally, a quantum system subjected to adiabatic 
driving acquires a geometric phase (or Berry phase) as 
well as a dynamic phase. Unlike the dynamic phase, the 
geometric phase depends solely on the trajectory of the 
parameters in the Hamiltonian, and thus is stable against 
certain types of fluctuations, which have been experimen- 
tally observed in various systems [ft-ill,]. However, the 
connection of the phase to the Landau-Zener-Stiickelberg 
(LZS) interferometry has not yet been experimentally 
demonstrated ,1^. 

The interferometry is composed of successive Landau- 
Zener (LZ) transitions, which occur at an avoided cross- 
ing of a driven quantum two-level system. The LZS inter- 
ferometry is similar to the Mach-Zehnder interferometry 
in optical systems in that the role of the LZ transition 
corresponds to that of a coherent beam splitter. The in- 
terference of the LZS interferometer is governed by the 
phase difference accumulated by the two instantaneous 
energy eigenstates between the subsequent transitions. 
The LZS interference has been observed in diverse phys- 
ical systems from atomic or optical systems 13l4l5[ to 
solid state systems |16l4l9| , where all of the observations 
are dominated by the dynamic phase or a combination 
of the dynamic and geometric phases [l3| . 

In our paper, we report the first experimental real- 



ization of the LZS interferometry controlled exclusively 
by the geometric phase inspired by the proposal of Ref. 



12l | in a single trapped ion system, which is capable of 



simulating other quantum two-level (qubit) systems. We 
observe the robustness of the geometric phase against 
immense operational errors in all possible control pa- 
rameters by artificially introducing noise into the sys- 
tem. Our demonstration of strong immunity sheds light 
on the possibility of examining the geometric phase in 
more complex systems which might be subject to large 
fiuctuations in control parameters. Furthermore, our re- 
alization contains the basic procedure for the adiabatic 
quantum simulation and can be extended to investigate 
and harness the geometric phase in many- body systems 
for quantum information processing |20l . |2l| . 

The geometric LZS interferometry of a qubit system 
can be described by the following Hamiltonian, 

i^GLZ it) = ^T • Beff (t) , (1) 

where a = ((Ta;,crj^,cr^) is the vector of Pauli matrices 
and the effective magnetic field Boff = {B^, By, B^) — 
{ilcosip, fl sin ip, A). Here A is the energy splitting of 
the qubit system at the avoided crossing and fl is the 
driving strength of the system. We denote H = 1 for 
convenience. 

We realize the geometric LZS interferometer in a sin- 
gle ^^^Yb^ ion as a model qubit. The single ^^-'^Yb^ ion 



is trapped in a four-rods radio frequency trap [22, |2 



We map the two internal levels of the ' Yb ion in 
the 5*1/2 ground state manifold to the qubit states, 
which is represented by: \F = l,mp — Q) = jt) and 
|F = 0, mp = 0) = 14,). The energy splitting of the two 
levels results from the hyperfine interaction and the tran- 
sition frequency between the |t) state and the jj,) state is 
whf = (27r) 12642.821 MHz. The coherent driving is im- 
plemented by the microwave with a frequency detuned by 
A(i) = ojhf — WM(i) and strength $7. We control A and 
r2 by mixing the microwave signal with the output of an 
arbitrary wave form generator with 1 GS/s, which is sig- 
nificantly high enough to ignore the sampling effect com- 
pared to typical operation time, a few hundred /iS. By 
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FIG. 1: (a) The control sequences of A and fi in the Bcff = (Jlcos;^, fisin;^, A) for the reahzation of the geometric LZS 
interferometry. The sequences are composed of: 1 preparation [ti,t2], 2 the 1°' LZ transition [f2,i3], 3a, 3b adiabatic evolutions 
for aquiring geometric phase [^3,^4] & [^6,^7], 4 the 2°'' LZ transition [fTjtg], 5 final rotation [igjfg] and spin echo sequence 
[t4,/;6l- At t = fe, we change the phase of <^ from to 0o- Adiabatic procedures are noted as red color (l,3a,3b,5). The colors 
and numbers used in (a) remain consistent in all the other figures, (b-d) The trajectories of the Beff following the control 
sequences, (e-g) The description of the geometric LZS interferometry in iJ — B^ space, where the hyperbolic curves indicate 
the adiabatic eigenenergies. In contrast to the standard energy diagram for the LZ transition, the geometric LZS interferometer 
should be described in 3-D space due to the phase information. The volume of the orange spheres corresponds to the population 
of the adiabatic eigenstates. (h) The trajectories of the Bcff for the adiabatic spin echo sequences, (i) The description of the 
spin echo sequences E — B^ig space. 



going to interaction picture defined by ifo = ^<^m(Oi ^^ 
obtain the geometric LZS Hamiltonian, i?GLZ (HJ , where 
Lp is the phase of the microwave source. In experiment, we 
first apply Doppler cooling and initialize the state to the 
14,) state by the standard optical pumping technique with 
99.1% efficiency 22,|23|. At the end of the experimental 
sequence, we measure the population of the \X) state by 
applying the fiuorescent detection scheme [22|, |23| . 

Fig. [TJa) shows the sequences of experimental con- 
trols in A and fi for the geometric LZS interferome- 
try. Fig. [IJb-e) show the trajectories of the Beg and 
Fig. mf-i) show the evolutions of qubit state in £■ — B^^ 
space according to changes of parameters in the Hamil- 
tonian, where the hyperbolic curves indicate the adia- 
batic eigenenergies E± = ±^Aq -|- Vfi. The sequences 
are composed of the five main procedures and adibatic 
spin echo: 1 [^1,^2] adiabatic prepartion of the instanta- 
neous ground state of the initial Hamiltonian, 2 [^2,^3] 
LZ transtion, 3a, 3b [^3, t4]&[t6,i7] adiabatic evolutions 
to accumulate geometric phase, 4 [^7,^3] LZ transtion, 
5 [tsj^g] adiabatic rotation to transfer the final state on 
the measurement basis Uz, and S1,S2 [^4,^5] spin echo 
sequence to eliminate the dynamic phase and the Stokes 
phase. The interference pattern of the geometric phase 
is observed by measuring the population of the upper 
eigenstate in the measurement basis. 

For the implementation of the geometric LZS interfer- 
ometer, the adiabatic evolution, where a parameter of a 
Hamiltonian slowly changes and the system follows the 



ground state of the instantaneous Hamiltonian, plays cru- 
cial role as the main methods for state preparation, state 
detection, state inversion and geometric phase genera- 
tion. Therefore, we carefully investigate the validity of 
adiabaticity in our experimental realization. We also per- 
form the experimental study of the LZ transition, where 
the separation probability between two instantaneous en- 
ergy eigenstates is controlled by the sweeping rate in the 
vicinity of the avoided crossing. Based on the experimen- 
tal confirmation of the validity of adiabatic evolution and 
of ability to control LZ transition, we perform the geo- 
metric LZS interferometry. 

Adiabatic Evolution The adiabatic evolution is used 
for the sequences of l,3a,3b,5 as well as S1,S2. We 
apply adiabatic evolution to rotate the Bcff field about 
an axis on the xy-plane by changing the amplitude 51 and 
the detuning A of the microwave in the following manner, 

n (<) = IBcff I sine (t) , A (t) - IBeffI cosB{t) , (2) 

where |Boff| is the magnitude of the effective magnetic 
field, 9 is linearly increasing in time (t) = {Of /To) t. 
We chose the changing rate Of /Ta — 7r/(200^s) for the 
|Beff| = (27r)50 kHz, which is small enough to satisfy 
the adiabaticity. With the rate, the initial ground state 
II) evolves as a blue curve shown on the Bloch sphere of 
Fig.[2l^a), where (p = 0. The time-dependency of the pop- 
ulations of the It) state are measured and compared to 
the numerical calculations [Fig. [2l^b)]. The difference in 
population between the experimental data and the ideal 
adiabatic evolution is no more than 5.5%. 
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FIG. 2: (a) The adiabatic evolution of the state. The traces of 
the system depending on the speeds of changing parameters 
are shown on the Bloch sphere in case of <^ = 0. The red 
curve is the trace of the perfect adiabatic evolution from the 
I4-) state to the ground state of J/glz = ^'^^ + '2'^^- We 
change A and fi as shown in Eq.((2]). The blue and the green 
curves are the trajectories with speeds Of /Ta = 7r/(200/is) 
and 7r/(60/is), respectively, with |Bcfi| = (27r)50 kHz. (b) The 
population of the W) state. During the adiabatic rotation, the 
rotating speed is 7r/(200^s), which is used though out the rest 
of the paper. The solid line shows numerical estimation of 
the evolution in W) state population and the + symbols show 
the experimental data after averaging the measurements of 
1500 trials, where the projection error of each point (0.012) 
is smaller than the size of the symbol. Note that we use the 
same coventions and perform the same number of trials for 
one data point throughout this paper. 



LZ Transition The LZ transition is used for the se- 
quences of 2 and 4. The LZ transition is investiaged 
including the time-resolved measurement of the tunnel- 
ing dynamics similarly to the demonstration with cold 
atoms [2J|. LZ tunneling occurs in the vicinity of the 
avoided crossing and in in the long time limit, the prob- 
ability transferred to the upper energy eigenstate of the 
adiabatic basis after the transition is characterized by 



Pi 



LZ 



exp 



2\v\ 



(3) 



where v = ^\ 



V = ^ Q_r,- In the adiabatic impulse approach, 
the transition in the adiabatic basis HV'-) j IV'+)} is de- 
scribed by the evolution matrix J/lz ~ UaNUb with 



N ^ 



ol-VS 



VI - Plz 



(4) 



where ips is the Stokes phase j25| and f/a(&) — 
exp (ifa^bjcr^), where ^^f^) is the dynamic phase accu- 
mulated after (or before) the LZ transition point. 

Experimentally, we use the sequences 1, 2 and 5 of 
Fig. [1] to study the LZ transtions: [Sequence 1] We 
prepare the ground state of the Hamiltonian ([T]) with 
fii = (27r)49.24 kHz, Aq = (27r)8.68 kHz and (p = 
by adiabatically rotating the {],) state about y-axis; [Se- 
quence 2] We change n{t) = (l — 2y) fJ; linearly in 
time; [Sequence 5] At time t, we adiabatically bring 



FIG. 3: (a) LZ tunneling probability Plz as the function of 
the total sweeping time T in the linear changes of Q, where 
Ao = (27r)8.68 kHz and fii = (27r)49.24 kHz. (b) The LZ 
transition dynamics for the three exemplary cases: the total 
durations T are 90 fis, 157 /is, and 300 fis, which provide 
0.5, 0.3 and 0.1 tunneling probabilities after the transitions. 
The oscillatory behaviors near the transition points [19l . |24| 
are clearly observed and precisely agreed with the numerical 
calculations. 



the state to that in the measurement basis {\ip- (t)) — )■ 
IJ,) , |i/)+ (t)) — >■ It)), which enables us to measure the 
population of the excited state |V'+ (t)) by observing the 
probability of the |t) state, which is equivalent to the 
transition probability. The population of the |t) state 
after the transition versus the duration of the transition 
is plotted in Fig. |31[a), where experimental results and 
the transition formula of Eq. (jH]) are in precise agree- 
ment. We also observe the transient dynamics and oscil- 
latory behavior of LZ tunneling in the vicinity of tran- 
sition point with various speeds of changing J7 shown in 
Fig. ^h) with improved quantitative agreements com- 
pared to the Ref. 19|, |24|. Note that there is no fitting 
parameter at the theoretical expectations in Fig. [31 since 
the parameters are independently measured. 

Geometric LZS Interferometer The sequences for 
the geometric LZS interferometer is described in Fig. 
[Tfa), where we use the squence 1,2 and 5 in the same way 
to the study for the LZ transition as described above. Af- 
ter the squence 1,2, the geometric phase is accumulated 
by the adiabatic evolutions shown in sequence 3a and 3b 
combined with the adiabatic spin echo sequence SI and 
S2 in Fig. [11 

[Sequence 3a,3b] & [Sequence S1,S2] The accu- 
lumated geometric phase of each adiabatic eigenstate 
{IV'-) , \i'+)} is given by. 



7± 



t3 



(V'± {t)\ ^ |V'± {t))dt 



(5) 



The difference of geometric phase between upper and 
lower eigenstates after the evolution is only determined 
by the rotating angle (po = (7+ — 7-) of the Beff shown 
in Fig. [ijc), which is independent of the energy differ- 
ence, the field amplitute change or the duration of the 
interferometry. The dynamic phase is also acquired in 
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FIG. 4: (a) The interference pattern of the geometric LZS interferometery. The final W) state populations are only determined 
by the geometric phase of the rotation angle ipQ. The solid black line comes from the theoretical expectation Eq. (|8]) and '+' 
symbols are experimental results with the average of 1500 trials, (b) The experimental demonstration and (c) the theoretical 
estimation of the immunities in the interference contrast against the errors in the amplitude Q, and the frequency A of Bcff . 
The amplitute and the frequency errors are scaled by the ideal Q. and the minimum frequency gap Aq. The experimental figure 
was constructed by the fringe contrasts of 1139 random choices of the amplitute Jl' with 1500 repetitions per choice. In our 
scheme, the constasts of the geometric LZS interferometer are clearly revealed even close to 30% and 50% fiuctuations of the 
frequency and amplitude. 



the adiabatic processes. In order to erase the dynamical 
phase, we adiabatically invert only the state as depicted 
in Fig. [TJeji) in the middle of the evolution [26|. Note 
that the Boff is conserved (same before and after the spin 
echo sequences). The time evolution operation Ug for 
these adiabatic stages including the spin echo sequences 
is written as, 



Ug = 



j«7- 





0*7+ 



(6) 



[Sequence 4] The second LZ transition at [^7,^3] in 
Fig. [Tfa) should be described in the basis of the inverted 
Beff, which results in (/lz = CTa^C/LZf"^. The final state 
\4>{tf)) after the second LZ transition can be expressed 
in the following equation. 



IV'/) = Ui^zUgUi^z \i) , 



(7) 



where the initial state is prepared at the instantaneous 
eigenstate with lower energy of the beginning Hamilto- 
nian. Note that the Stokes phases occurring at the first 
and second LZ transitions are effectively cancelled out 
because of the inversion of the state. Therefore, the final 
result is insensitive to the fluctuation of the Stokes phase 
as well as the dynamic phase. 

[Sequence 5] At the final rotation [isj ^9] in Fig. [IJa), 
we adiabatically transfer the population of the adiabatic 
basis to the measurement basis as discussed in the adi- 
abatic process. Finally the population of the |t) state 
after the geometric LZS interferometry can be described 
by the simple formula, 



P\^) = 2Plz (1 - Plz) (1 - cos(^o) . 



(8) 



In the experiment, we set the transition probability 
-Plz = 1/2, which simplifies ([S]) to P\-^^ = 1 (1 — co&Lpo). 



Fig. Ufa) clearly shows that the populations of the |t) 
at the end of the interferometry are solely determined by 
the geometric phase acquired during the adiabatic evo- 
lution, which is exactly the rotating angle <^o- Fig- Sfb) 
and (c) show the immunity of our interferometry to the 
errors in the amplitude J7 and the frequency A that are all 
parameters in the effective field Boff. For the amplitude 
error scaled by the ideal strength, about ±50% changes 
reduce the contrast of the interference pattern by 20%, 
and for the frequency errors relative to the minimum gap 
Aq, ±30% offsets decrease the contrast by 20%. Here, we 
assume 11' and A + (5 to be unchanged in a single experi- 
ment. In this case, the dynamic phase and Stokes phase 
are always cancelled out, while the probability of the LZ 
transition changes between different experiments, which 
cause the reduction of the interference contrast. We note 
it has been shown that the geometric phase is robust 
against the fast fluctuations of the control parameters, 
where the evolution time is longer than the typical noise 
correlation time 



Hi 



In conclusion, we have realized a clear connection be- 
tween the geometric phase and the LZS interferometry 
and have observed the interference of pure geometric 
phase. We have demonstrated the robustness of the in- 
teference against variations of all the parameters in Bcfj , 
which shows the possibility of observing such interef- 
erence even in more complex systems, including super- 



conducting qubit systems [12|, |26[ , quantum dot systems 



27| , NV center diamond systems [19|, etc. Within 



trapped ion system, our research makes possible a new 
direction in performing multi-qubit operations with the 
geometric phase in internal states [20|, |28| , which can pro- 
vide strong robustness and high fidelity of operation be- 
yond the current methods that use the geometric phase 
in harmonic oscillators 



29l 30f . We also note that the ex- 



perimental method used in the realization can be directly 



applied to adiabatic quantum computation. 
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